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1 INTRODUCTION 


ABSTRACT 

The third data release of Gaia has provided low resolution spectra for ~ 100 000 white dwarfs 
(WDs) that, together with the excellent photometry and astrometry, represent an unrivalled 
benchmark for the study of this population. In this work, we first built a highly-complete 
volume-limited sample consisting in 12718 WDs within 100 pc from the Sun. The use of 
VOSA tool allowed us to perform an automated fitting of their spectral energy distributions 
to different atmospheric models. In particular, the use of spectrally derived J-PAS photometry 
from Gaia spectra led to the classification of DA and non-DA WDs with an accuracy > 90%, 
tested in already spectroscopically labelled objects. The excellent performance achieved was 
extended to practically the whole population of WDs with effective temperatures above 5500 K. 
Our results show that, while the A branch of the Gaia WD Hertzsprung-Russell diagram is 
practically populated by DA WDs, the B branch is largely formed by non-DAs (65%). The 
remaining 35% of DAs within the B branch implies a second peak at ~ 0.8 Mọ in the DA- 
mass distribution. Additionally, the Q branch and its extension to lower temperatures can be 
observed for both DA and non-DA objects due to core crystallisation. Finally, we derived a 
detailed spectral evolution function, which confirms a slow increase of the fraction of non-DAs 
as the effective temperature decreases down to 10500 K, where it reaches a maximum of 36% 
and then decreases for lower temperatures down to ~ 31%. 


Key words: white dwarfs — stars: evolution — Galaxy: stellar content- astronomical data bases: 
miscellaneous — catalogues — virtual observatory tools 


gravity by Fermi electron pressure and they are doomed to a long- 
standing cooling process. As this physical process is relatively well 


The excellent quality of astrometric and photometric data provided 
by the ESA mission Gaia has been recently improved by the publica- 
tion in its third data release (DR3) of nearly 200 million spectra (De 
Angeli et al. 2022). These spectra, although having a very low res- 
olution (R ~ 60), represent an invaluable source of information for 
a wide range of stellar objects in our Galaxy. In particular, nearly 
100000 spectra correspond to white dwarfs (Gaia Collaboration 
et al. 2022a). 

As it is well known, white dwarfs are stellar remnants of low- 
to intermediate-mass main sequence stars (Althaus et al. 2010). 
Since nuclear fusion processes have ceased after the earlier stages 
of their lives, white dwarfs are stellar objects supported against 


* E-mail: fran.jimenez-esteban @cab.inta-csic.es 


© The Authors 


understood and some white dwarfs can reach very old ages (10 Gyr 
or older), they can be used as reliable cosmochronometers (Fontaine 
et al. 2001), revealing valuable information about the history and 
evolution of our Galaxy (Garcia-Berro & Oswalt 2016). 

The source of energy in the deep interior of white dwarfs — 
basically due to the gravothermal energy released by the ions and 
eventually provided by core crystallisation, phase separation, and 
sedimentation of minor species among other processes (see Isern 
et al. 2022, for a recent review) — is controlled by a thin partially 
degenerate layer through which the heat is radiated away. In the 
canonical model, this layer is formed by helium with a mass around 
107? Mo, representing less than 2% of the total white dwarf mass, 
and in most of the cases (~ 80%), an extra thinner layer of hydrogen 
with a mass between 10715-1074 Mo that lies on top of the helium 
one. 
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Observationally, characterisation of white dwarfs has been car- 
ried out mainly from spectroscopy, based on the atmospheric ob- 
served features, for instance through the identification of Balmer and 
helium lines, among others (Sion et al. 1983). If the white dwarf 
spectrum presents hydrogen lines the white dwarf is labelled as 
DA, while if the spectrum presents absorption helium lines, He1 or 
He 1, the white dwarf is named DB or DO, respectively. Although 
less common, it is possible to identify some metals within white 
dwarf atmospheres, such as carbon, named DQ, or other heavy ele- 
ments such as Ca n or Fe n, assigned to the DZ spectral type. Some 
white dwarfs can also present featureless spectra, named in this 
case as DC. This last condition is usually satisfied by cool objects, 
in particular for hydrogen-dominated atmospheres with effective 
temperatures (Tef) below ~ 5000 K, or for helium-dominated atmo- 
spheres with Teg £ 11 000K. This initial list of spectral types has 
been extended to include mixed cases. Thus, we can find hydrogen- 
dominated atmospheres with some helium features, named DAB, or 
the contrary case, DBA, or any possible combination such as DBQ, 
DAZ, and so on. Also possible is the detection of a magnetic field 
or the existence of some variability in the white dwarf spectrum, 
adding in these cases and extra H or V, respectively, to the basic 
spectral type, e.g., DAH, DAV, DQV. All in all, we can summarize 
the spectral classification of the white dwarf population in two main 
classes: those with a hydrogen-dominated atmosphere (which is by 
far the most common case), DA, and the rest, named as non-DA. 
In any case, the importance of determining the atmospheric content 
of a white dwarf is crucial, not only for understanding the physical 
processes in their evolution, but also to derive reliable individual 
parameters such as masses, temperatures, and ages. 

Spectral identification has historically relied on visual inspec- 
tion of spectra. However, with the advent of large spectroscopic sur- 
veys, such as Sloan Digital Sky Survey (SDSS; York et al. 2000), 
Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAM- 
OST; Zhao et al. 2012), Gaia (Gaia Collaboration et al. 2016) and 
many others, new automated data analysis techniques are required. 
In this sense, the Spanish Virtual Observatory (SVO) developed a 
powerful tool, the Virtual Observatory Spectral energy distribution 
Analyzer (VOSA!; Bayo et al. 2008). VOSA is a Virtual Observatory 
(VO) tool that allows the user to build spectral energy distributions 
(SEDs) from both private and public photometric data within the 
VO, of thousands of objects at a time, and to derive their physical 
properties by comparing the observed SEDs to different collections 
of theoretical models. Moreover, it can estimate the goodness of the 
fitting of the model to the observed data. 

In this work we aimed at thoroughly analyzing the 100 pc 
nearly-complete volume-limited sample of white dwarfs with the 
main objective of individually classifying them based on their main 
spectral type class (DA or non-DA). The paper is, thus, organized 
as follows. After this introduction, in Section 2 we describe our 
methodology for building our volume-complete observed sample 
of white dwarfs within 100 pc from the Sun. In Section 3, we detail 
the different white dwarf atmospheric models used, and how the 
SEDs were built from VO photometry or from synthetic photometry 
derived from the Gaia spectra. In Section 4, we introduce a set of 
white dwarf spectral type estimators based on the VOSA fitting 
and also, for comparative purposes, others estimators based on data 
found in the literature. In Section 5 we apply the best estimator found 
on the basis of an already classified population of white dwarfs to 
the entire population, thus obtaining reliable distributions of stellar 
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parameters, such as mass and effective temperature. Finally, we 
summarize our main results and conclusions in the last section. 


2 THE GAIA-DR3 100 PC WHITE DWARF SAMPLE 


In this section we describe our methodology for identifying white 
dwarf candidates within 100 pc in the Gaia-DR3 catalogue. 

In Jiménez-Esteban et al. (2018) (hereafter JE18), we demon- 
strated that the completeness of the 100 pc white dwarf sample 
selected using the second Gaia data release (DR2; Gaia Collabora- 
tion et al. 2018) reached ~ 94% for a parallax relative error lower 
than ~ 10%. We also noted that smaller samples in volume are not 
necessarily more complete due to the Lutz-Kelker bias (Lutz & 
Kelker 1973) and that for larger volumes, the completeness also 
decreases due to the magnitude limit of the Gaia satellite (Gaia 
Collaboration et al. 2016). As a consequence, we focused in this 
work on the construction of a white dwarf catalogue up to 100 pc 
by selecting Gaia-DR3 sources with parallaxes larger than 10 mas, 
considering the errors, and imposing relative errors smaller than 
10% in parallax and in BP- and RP-band photometry. Photomet- 
ric errors in the BP and RP bands were especially important for 
faint sources in crowded areas. Since the G-band photometry had 
lower error than the other two bands, we did not impose any cut in 
the G-band. Thus, we queried the Gaia-DR3 catalogue” using the 
following criteria: 


© w-30y = 10mas and w/a, = 10 
e FBP/CTFpp = 10 and FRP/C Fpp = 10 
e RUWE<1.4 


The cut in the RUWE (Renormalised Unit Weight Error) pa- 
rameter prevented against poor astrometric solutions (Lindegren 
et al. 2021). These criteria guaranteed that the sources were indeed 
within the 100 pc volume. In addition, we estimated the corrected 
BP and RP flux excess factor (C*) and its scatter (7+), and applied 
a 30c» cut following the recommendations by Riello et al. (2021). 
This criterion excluded some objects in the coolest end of the white 
dwarf sequence, where the oldest white dwarfs are expected. How- 
ever, this guaranteed the consistency between the G-band and the 
BP and RP photometry. 

We show in Fig. 1 the Gaia Hertzsprung—Russell (HR) diagram 
for all sources within 100 pc with good photometric and astrometric 
data, as defined by the above criteria. Since the BP-band photom- 
etry has a known bias for red and faint sources (Fabricius et al. 
2021), we used the G — Grp colour instead of the more commonly 
adopted Gpp — Grp. The 10% of uncertainty of the Gaia fluxes 
corresponds to an uncertainty of ~ 0.14 mag in the Gaia colour. On 
the other hand, the 10% uncertainty in the parallax corresponds to 
an uncertainty of ~ 0.2 mag in the absolute G-band magnitude. This 
diagram clearly shows the Main Sequence, the Red Giant Branch, 
and the White Dwarf Sequence. Based on the location within the 
colour-magnitude diagram, we selected 12718 Gaia-DR3 sources 
located in the White Dwarf Sequence (blue and red dots in Fig. 1), 
thus representing our 100 pc sample. It is worth noting that, 7480 
of these sources (blue) were already in the Gaia-DR2 white dwarf 
catalogue? by JE18, and that 5238 (red) are new incorporations 
from the Gaia-DR3 catalogue. 

The Gaia-DR2 white dwarfs are located in Fig. 1 at bluer colour 
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Figure 1. Gaia Hertzsprung—Russell diagram of all Gaia sources within 
100 pc considered in this work. White dwarfs already identified in the Gaia- 
DR2 catalogue by Jiménez-Esteban et al. (2018) are plotted in blue, and 
in red the new candidates identified in the Gaia-DR3 catalogue. Those 
stars fulfilling our 100 pc selection criteria but not selected as white dwarf 
candidates are shown in grey. 


(G—Gprp <0.5 mag), while most of the newly identified have redder 
colours up to G — Grp ~1.2 mag. This is because JE18 imposed 
a colour limit of Ggp — Grp < 0.8 mag, corresponding to an 
Te of ~ 6000 K, to avoid contamination by field stars with bad 
Gaia-DR2 astrometry. Thanks to the improvement in the astrometry 
measurements of the Gaia-DR3, this colour cut is no longer needed, 
allowing us to study the red part of the White Dwarf Sequence. There 
is also a number of hotter white dwarfs (Ggp— Grp < 0.8 mag) that 
were not included in the Gaia-DR2 white dwarf Catalogue by JE18. 
51 of them are located in the HR diagram above the box defined in 
that work. The rest was not included probably due to lower quality 
of the Gaia-DR2 astrometry and/or photometry. On the other hand, 
there are 139 white dwarf candidates within 100 pc identified by 
JE18 from the Gaia-DR2 Catalogue which are now excluded from 
the new catalogue using Gaia-DR3. This is probably due to a wrong 
astrometric solution in the Gaia-DR2 catalogue. 

In conclusion, the catalogue of white dwarfs within 100 pc pre- 
sented in this work, not only extends to redder colours than the pre- 
vious one, but also discards misidentified objects and includes blue 
missing sources. This catalogue, contining the basic information of 
the selected sources, including the spectral classification (Sect. 4) 
and the physical parameters (Sect. 5.2) derived in this work, can 
be accesed electronicaly through our online catalogue service (see 
Appendix A). 


2.1 Contamination 


We used the SIMBAD database to evaluate the contamination in 
our white dwarf catalogue by other types of sources. We searched 
for SIMBAD counterparts within 3” of our Gaia white dwarf can- 
didates using the Gaia-DR3 coordinates translated to epoch J2000 
and corrected from Gaia-DR3 proper motions. Approximately 98% 
(12 454) of the sample had a counterpart in SIMBAD, all of them 
with a classification, at least tentative. The results are shown in 
Table 1. If we assume the SIMBAD classifications to be correct, 
including those flagged as candidates, and taken into account the 
position of these sources in the Gaia colour-magnitude diagram 
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Table 1. Available SIMBAD classification for 12 454 of our selected Gaia- 
DR3 white dwarfs. Classifications in agreement with the white dwarf nature 
of our selected candidates are listed above and shown in italics. The estimated 
contamination is < 0.1%. 


SIMBAD class N (Candidates) 

white dwarf 3640 (7242) 

Blue object 4 

Star 80 

High proper motion 1449 

In Binary / Multiple system 18 
Hot subdwarfs 12 (2) 
Cataclysmic Variable 4 (1) 
Variable star 1 
Emission-line star 1 


(see Fig. 1), one can conclude that most (>99.8%) of them have 
classifications compatible with white dwarf nature. 

The highest source of contamination seems to be hot subdwarfs 
(~ 0.1%) and cataclysmic variables (~ 0.04%). However, although 
hot subdwarfs have similar colours than hot white dwarfs, they 
are more luminous and should occupy a different region in the 
HR diagram. So, the classification provided by SIMBAD is likely 
to be incorrect, and these objects are most probably misclassified 
white dwarfs. Thus, the expected contamination in our white dwarf 
catalogue is very low (< 0.1%). 


2.2 Completeness 


Our 100 pe white dwarf sample is volume-limited and has a very 
low degree of contamination sample. However, it is subject to the 
astrometric and photometric constrains that were imposed to the 
Gaia-DR3 catalogue as explained above. Recently, from the Gaia- 
DR2 catalogue Hollands et al. (2018) estimated a space-density of 
4.49 x 1073 pe~3 for the white dwarf population within 20 pe from 
the Sun. Assuming the same space-density upto 100 pc, we would 
expect ~18.800 white dwarf in this volume. This implies that the 
completeness of our sample is ~70% of the entire 100 pc white 
dwarf population. This estimation is in agreement with Torres et al. 
(2022). 

Nevertheless, the level of completeness is not homogeneous 
along the Gaia colour range presented by the sample. Fig. 5 of 
Jiménez-Esteban et al. (2018) shows the completeness as a function 
of the Ggp — Grp colour of a similar sample within 100 pc obtained 
from the previous Gaia-DR2 catalogue. The completeness is almost 
100% for Ggp — Grp <Omag and it continuously decreases to 
~70% for the reddest colours Ggp — Grp = 0.8 mag. It is expected 
that the completeness continues decreasing for redder colours than 
Gpp — Grp = 0.8 mag. 

Our sample was built from the Gaia-DR3 catalogue, which has 
both better astrometry and photometry than the previous Gaia-DR2, 
so the completeness of the present sample is expected to be higher 
than the previous one. Consequently, we conclude that most of the 
missing sources are at the reddest colours, and the completeness of 
the spectroscopically classified sources (Ggp — Grp < 0.86 mag; 
see Sect.4) is larger than 90%. 

Regarding the missing objects, one cause of missing sources is 
the confusion in the Galactic plane. Using our white dwarf sample, 
we compared the number of white dwarfs per square degree in the 
Galactic plane (|b| < 10 deg) and out of it. We found that the white 
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dwarf sky density in the Galactic plane is 5% lower than in the rest 
of the sky. Thus, we estimated ~120 missing white dwarfs due to 
confusion. This is ~1% of our sample. 

Another cause of missing white dwarfs is the number of double 
degenerate systems (DWD) which were not resolved by Gaia and 
were counted as single white dwarfs. Torres et al. (2022) used 
different models to estimate the fraction of such systems. The best 
model estimated a fraction between 1% and 3% (depending on the 
common-envelope treatment) of unresolved DWDs within 100 pc. 

Finally, a third cause of missing sources are the binary systems 
made up by a white dwarf and a main sequence AFGK star, also 
known as Sirius-like systems. Holberg et al. (2013) estimated that 
8% of the white dwarfs in the solar vecinity are members of such 
systems. Since many of them are now resolved by Gaia, we can 
used this value as an upper limit of the missing white dwarfs in our 
sample. 

Summing up the contribution of missing sources, we conclude 
that less than 10% of white dwarfs located within 100 pc are not part 
of our catalogue in the colour range Ggp—Grp < 0.86 mag, while in 
the reddest end of the white dwarf sequence completeness decreases 
and the missing source would be mainly due to the limitations of 
Gaia capabilities in both astrometry and photometry. 


3 SPECTRAL ENERGY DISTRIBUTION ANALYSIS 


In this section we present the analyse of the SEDs of our 100 pc 
white dwarf sample built from public photometric data and from 
Gaia spectra. Our purpose was to classify them in two main spec- 
tral classes: DA and non-DA. To that end, we used two different 
grids of white dwarf atmosphere models: hydrogen-dominated at- 
mospheres (type DA) and helium-dominated atmospheres (note that 
even though DB models were used, others spectral types were also 
expected and we therefore named this group as type non-DA). We 
did not account for interstellar extinction since our sample is located 
at short distances. 


3.1 White dwarf atmospheric models 


The collection of models used in the analysis adopt LTE (local ther- 
modynamic equilibrium), hydrostatic equilibrium, plane-parallel, 
one-dimensional structure, and convection with the mixing-length 
approximation where appropriate. The DA models consider pure hy- 
drogen composition, and the non-DA models consider helium with 
a small trace of hydrogen (log N(H)/N(He) = -6). The version of the 
mixing length is ML2 (Tassoul et al. 1990), with the ratio of mixing- 
length to pressure scale height 0.7 in the DA and 1.25 in the non- 
DA white dwarfs. Basic methods and data are described in Koester 
(2010). More recently many improvements were implemented, the 
most important being: improved treatment of molecules, non-ideal 
effects, unified line broadening theories for the strong lines of Mgl, 
Mgl, Cal, Call, and Lyman alpha. Also, the hydrogen Stark profiles 
by Tremblay & Bergeron (2009) and Tremblay (2015, priv. comm) 
are included. 

The DA grid contains 1260 model spectra and covers effective 
temperatures from 3000 to 20 000 K in steps of 250 K (note that be- 
low ~ 5000 K these white dwarfs are DCs), from 20 000 to 30 000 K 
in steps of 1000 K, and from 30 000 to 40 000 K in steps of 2 000 K. 
For each effective temperature, the surface gravities range from 6 to 
9.25 dex in steps of 0.25 dex. The non-DA grid contains 666 model 
spectra and, in this case, the effective temperatures cover the 5500 
to 20000 K range in steps of 250 K (below ~ 10000 K these white 


dwarfs are also featureless DCs), the 20000 to 30000 K range in 
steps of 1000 K, and the range from 30000 to 40 000 K in steps of 
2 000 K. For this grid, the surface gravities range from 7 to 9 dex in 
steps of 0.25 dex for each effective temperature value. 


3.2 SED from VO archives 


To build the SEDs of the 100 pc white dwarf sample from public 
archives within the VO, we collected data from different photomet- 
ric catalogues, which are listed in Table 2. Detailed information of 
the adopted filters can be found using the SVO Filter Profile Ser- 
vicet (Rodrigo et al. 2012; Rodrigo & Solano 2020). In order to 
avoid misidentifications, the Gaia coordinates for all sources were 
corrected from proper motion and transformed to an epoch closer to 
the observation epoch for each of the queried survey. A 3” search 
radius was adopted. This allowed us to build the observational SEDs 
from the UV to the mid-IR wavelength range for most of the sources. 

For 463 white dwarf candidates, no additional photometry to 
Gaia data were found, so their SEDs could not be fitted with VOSA. 
For the other 12 255 sources, each individual SED was fitted to both 
DA and non-DA white dwarf synthetic spectra outlined in Sect. 3.1. 
Of them, we obtained a reliable fit, defined as Vgfp < 155 , for 10 447 
sources when using either DA or non-DA model spectra. 1104 white 
dwarfs obtained good fits only when using DA model spectra, while 
65 white dwarfs obtained good fits only when using non-DA model 
spectra. In total, we obtained a good fit for 11 616 sources (~ 95% 
of the fitted sample). 

For illustrative purposes, in the left panels of Fig.2 we show 
the observational SEDs built from VO photometric catalogues (red 
circles), together with the synthetic photometry (blue squares) and 
the synthetic spectra (grey solid line) that best fit the observations. 
The two upper rows correspond to two white dwarfs which best-fit 
model is a DA, and the two bottom rows to two white dwarfs which 
best-fit model is a non-DA, with different effective temperatures and 
surface gravities. In the right panels of Fig. 2 we also show for the 
same sources the SEDs built with the J-PAS photometry obtained 
from Gaia spectra (see Sect. 3.3). 


3.3 SED from Gaia spectra 


The new Gaia-DR3 has provided low-resolution spectra for more 
than 200 millions sources (De Angeli et al. 2022). The spectra cover 
the optical to near-infrared wavelength range, from 330 to 1050 nm, 
approximately. 12 342 (97% of the whole ample) white dwarfs in 
our sample had Gaia low-resolution spectra available. We used the 
Python package GaiaXPy® to construct the photometric SEDs of 
these sources from Gaia spectra using the Javalambre-Physics of 
the Accelerating Universe Astrophysical Survey (J-PAS; Benitez 
et al. 2014) filter system (Marin-Franch et al. 2012). This system is 
composed of 54 overlapping narrow-band (FWHM ~ 145 A) filters 
covering from 3780 to 9100 A, plus 2 broad-band filters at the blue 
and red end of the optical range, and 4 complementary SDSS-like 
filters. In total, the 60 J-PAS filters provide a low-resolution (R ~ 60) 


4 Available at http: //svo2.cab.inta-csic.es/theory/fps/. 

5 Vgfp: Modified reduced y”, calculated by forcing o-(Fops) to be larger 
than 0.1 Fobs, where o (F ops) is the error in the observed flux (Fops). This 
can be useful if the photometric errors of any of the catalogues used to build 
the SED are underestimated. Vgf; smaller than 10 - 15 is often perceived as 
a good fit. 

6 nttps://www.cosmos.esa. int/web/gaia/gaiaxpy 
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Figure 2. Spectral energy distributions of four distinct white dwarfs, each row presents an object (from the top to the bottom: Gaia-DR3 3750749378584 132992, 
Gaia-DR3 3020668542435696512, Gaia-DR3 2780434524599787136, and Gaia-DR3 5869567658943170048), illustrating DA and non-DA models with 
different effective temperatures and surface gravities. The red filled circles show the photometry available in VO archives (left panels) and the J-PAS 
photometry built from Gaia spectra considering all coefficients (right panels). The best-fit models are presented as grey solid line, where the blue squares are 
the synthetic photometry (calculated from the model) in each filter. The dashed line represents the Rayleigh-Jeans law adopted for longer wavelengths, when 


they are not covered by the models. 
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Table 2. Photometry data from VO archives used to build the SEDs. The detailed information on each photometric filter can be obtained from the SVO Filter 
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Profile Service (Rodrigo et al. 2012; Rodrigo & Solano 2020). 


Survey Facility/Instrument Spectral range Filters Reference 
GALEX UV GR6+7 GALEX UV FUN, NUV 1 
XMM-SUSS4.1 XMM/OM UV, visible UVW2, UVM2, UVW1, U, B, V 2 
Gaia DR3 Gaia visible Gpp, G, Grp 3 
APASS DR9 — visible B, v, g,r,i 4 
J-PLUS DR1 OAJ/T80Cam visible uJAVA, J0378, J0395, J0410, J0430, gSDSS 5 
J0515, rSDSS, J0660, iSDSS, J0861, zSDSS 

Tycho-2 TYCHO visible B, V 6 
SDSS DR12 SLOAN visible u, g,T,i, Z aT, 
IPHAS DR2 INT/WFC visible gR, Ha, gI 8 
VPHAS+ DR2 Paranal/OmegaCAM visible u, g, Ha,i,z 9 
Pan-STARRS1 DR2 PAN-STARRS visible, near-IR g.1r,i,Z,y 10 
Dark Energy Survey DR1 CTIO/DECam visible, near-IR g. r,i, z, Y 11 
DENIS — visible, near-IR I,J, Kg 12 
VISTA Paranal/VIRCAM visible, near-IR Z, Y,J,H, Ks 13 
UKIDSS UKIRT/WFCAM visible, near-IR Z, Y,J,H,K 14 
2MASS 2MASS near-IR J, H, Ks 15 
GLIMPSE I+1+3D Spitzer/IRAC mid-IR I1, I2, 13, 14 16 
AlIWISE WISE mid-IR W1, W2, W3, W4 17 


References: [1] Revised catalog of GALEX UV sources (GALEX UV GR6+7; Bianchi et al. 2017); [2] XMM-Newton Serendipitous Ultraviolet Source 
Survey catalogue (KMM-SUSS4.1; Page et al. 2012); [3] Gaia DR3 (Gaia Collaboration et al. 2022b); [4] AAVSO Photometric All-Sky Survey (APASS 
DR9; Henden et al. 2015); [5] Javalambre Photometric Local Universe Survey (J-PLUS DR1; Cenarro et al. 2019); [6] Tycho-2 Catalogue (Høg et al. 
2000); [7] Sloan Digital Sky Survey Photometric Catalogue (SDSS DR12 Alam et al. 2015); [8] INT/WFC Photometric H-Alpha Survey of the Northern 
Galactic Plane Catalogue (IPHAS DR2 Barentsen et al. 2014); [9] VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge (VPHAS+ 
DR2; Drew et al. 2016); [10] Panoramic Survey Telescope and Rapid Response System (Pan-STARRS1 DR2; Chambers et al. 2016); [11] Dark Energy 
Survey DR1 (Dark Energy Survey DR1; Dark Energy Survey Collaboration et al. 2016); [12] Deep Near Infrared Survey of the Southern Sky third release 
data (DENIS Denis 2005); [13] Visible and Infrared Survey Telescope for Astronomy (VISTA Cross et al. 2012); [14] UKIRT Infrared Deep Sky Survey 
(UKIDSS Hewett et al. 2006); [15] 2MASS All-Sky Catalog of Point Sources (2MASS; Skrutskie et al. 2006); [16] Galactic Legacy Infrared Midplane 
Survey Extraordinaire Source Catalog (Spitzer Science 2009, GLIMPSE I+II+3D); and [17] Wide-field Infrared Survey Explorer AllWISE Data Release 


(AIIWISE Cutri et al. 2021). 


spectrum (Bonoli et al. 2021) similar to Gaia. However, due to the 
low signal-to-noise ratio of the Gaia spectra at bluer wavelengths 
(Gaia Collaboration et al. 2022a), we did not use the 4 J-PAS filters 
with effective wavelength shorter than 4000 A, which reduces to 
56 the number of J-PAS filters available for building the SEDs. 
Detailed information on the J-PAS filter system can be found at the 
SVO Filter Profile Service’. 

For each source in our sample with a Gaia spectrum, we built 
two SEDs using the J-PAS photometric system, one taking into 
account all the coefficients of the spectrum (GJP) and one taking 
into account only the relevant coefficients as provided by Gaia- 
DR3 (GJP-trunc) (De Angeli et al. 2022). We imposed a threshold 
in the photometric error of 10% to each individual photometric 
measurement obtained with GaiaXPy. Thus, although most of the 
SEDs have 56 photometric points, in some noisy spectra the number 
of points is lower. For illustrative purposes, in Fig.3 we show two 
SEDs built considering all coefficients for a source with one of the 
highest and another with one of the lowest signal-to-noise ratio. 

Of the 12 342 Gaia spectra, we obtained 11 447 GJP SEDs with 
more than 4 photometric points, the minimum number necessary 
to be fitted by VOSA, and 12324 GJP-trunc SEDs. Of them, we 
obtained a reliable fit (Vgf, < 15) for 11 153 GJP and 11 639 GJP- 
trunc SEDs, when using either DA or non-DA model spectra. 290 
GJP and 678 GJP-trunc SEDs obtained good fits only when using 
DA model spectra, while no SED obtained a good fit only when 


7 nttp://svo2.cab.inta-csic.es/theory/fps/index. php? 
mode=browse&gname=0AJ&gname2=JPAS 


using non-DA model spectra. In total, we obtained a good fit for 
11443 GJP and 12317 GJP-trunc SEDs, almost all the analyzed 
SEDs obtained from Gaia spectra. To illustrate, we show some 
examples of the obtained GJP SEDs and their best fitting model in 
Fig. 2 (right panels). 


4 WHITE DWARF SPECTRAL CLASSIFICATION 


In this section a set of white dwarf spectral class estimators is 
built. Most of these estimators were based on the result of the 
fitting of the SEDs with certain white dwarf atmospheric models 
as described in the previous Section. The performance of these 
estimators was validated through the spectroscopically confirmed 
white dwarfs from the Montreal White Dwarf Database (MWDD; 
Dufour et al. 2017). 


4.1 Spectral estimators 


Following the VOSA procedure, in which each SED was fitted to 
two models (a DA and a non-DA model), and for which VOSA 


derived two reduced y?-values ( X5 a and KaD a> respectively), 


we built an estimator, i, that measures the probability of being a 
i 


hydrogen-dominated atmosphere (DA) white dwarf, Pha. as 
2 2 
pi = 1 | Xnon-DA ~ XDA +1 (1) 
DA 9 2 +y 
Xnon-DA * XDA 
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Figure 3. Spectral energy distributions of two white dwarfs. Gaia-DR3 
368351950388 1169920 (upper panel) has one of the highest signal-to-noise 
ratio and blue colour (Ggp — Grp ~ —0.32 mag). On the contrary, Gaia-DR3 
5174110233491949312 (lower panel) has one of the lowest signal-to-noise 
ratio and red colour (Ggp — Grp ~ 0.86 mag). The filled circles show the 
J-PAS photometry built from Gaia spectra considering all coefficients. In 
red, those photometric points with a relative flux error lower than 10% are 
shown, and then considered in the model fit, and in black are shown those 
not considered due to larger errors. The Gaia spectra are presented as grey 
solid lines. 


Hence, we classified an object as a DA white dwarf if ae A Zm 
where 77 is a threshold value, otherwise we considered it as a non- 
DA. In our analysis we adopted an initial value of 7 = 0.5. 

This estimator can be generalized by using any other y2-values 
available in the literature. Hence, up to 6 estimators were analyzed 
in this work. A description of them is presented as follows, while a 
summary is listed in Table 3. 


4.1.1 The VOSA, VOSA-GJP, and VOSA-GJP-trunc estimators 


In Section 3.2 we described the use of VOSA for analyzing white 
dwarf SEDs. For those objects with a reliable fit, i.e., Vgfb<15, us- 
ing both DA and non-DA atmospheric model, and by means of eq. 1, 
we built three different estimators. The first one, called VOSA, cor- 
responds to the case that the SEDs were built from public archives 
within the VO (see Section 3.2). In the case that the SEDs were 
derived from Gaia spectra (see Section 3.3), we defined two extra 
estimators: if all the coefficients of the white dwarf spectra were 
taken into account, our estimator is named VOSA-GJP; and in the 
case that the truncated version of the coefficients was adopted, our 
estimators is called VOSA-GJP-trunc. 
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4.1.2 The GF21-I and GF21-II estimators 


The analysis of the Gaia Early Data Release 3 white dwarf pop- 
ulation performed by Gentile Fusillo et al. (2021) also provides a 
fitting of different theoretical atmospheric models to the three Gaia 
passbands: G, Grp, and Gpp. Three models were used: one for 
hydrogen-pure atmospheres of Tremblay et al. (2011) with Lyman 
a opacity of Kowalski et al. (2010), another for pure-helium atmo- 
spheres based on Bergeron et al. (2011) models, and the last one for 
mixed H/He atmospheres with a ratio log N(H)/log N(He) = 1075 
based on models of Tremblay et al. (2014) and McCleery et al. 
(2020). All the three models cover the range of effective tempera- 
tures and surface gravities of the white dwarfs analyzed here (see 
Gentile Fusillo et al. 2021, and reference therein for further details). 

For each atmospheric model — hydrogen-pure, helium-pure, 
and mixed atmospheres — a y2-value, AM Kiepie Xfe-pure and 


j xe g respectively, was provided in Gentile Fusillo et al. (2021). 


From the first two y?-values, apie and Xfie-pure Gas and 
oe a> respectively), we build from eq. 1 our fourth DA prob- 
ability estimator called GF21-I. Similarly, from the first and third 
y2-values, Ne ies and XŽixe al YD a and Ne as a> respectively), 
we generated our fifth estimator named GF21-II. 


4.1.3. The MON22 estimator 


All except one of the estimators analyzed in this work follow eq. 1. 
The exception is the sixth estimator used in this work. It is rep- 
resented by the probability of being a DA estimated by Gaia Col- 
laboration et al. (2022a). In that case, that probability was directly 
derived from the application of a Random Forest algorithm to the 
J-PAS synthetic photometry derived from Gaia-DR3 spectra. 


4.2 Validating sample: the Montreal White Dwarf Database 


MWDD3 is an open access tool containing spectroscopically classi- 
fied white dwarfs published in the literature. At the time of writing 
this article the database contained 68 364 objects, 3098 of which are 
within 100 pc from the Sun. From them, we rejected those which 
are binaries, those with circumstellar disc, and those with tentative 
spectral clasification. This reduced the sample to 2886 sources. We 
used this data base to construct a validating sample. 

Itis worth saying here that our white dwarf atmospheric models 
cover a different range of temperatures, i.e., down to 3000 K for DA 
and down to 5500 K for non-DA (see Set. 3.1). To ensure a proper 
classification of the Gaia sample we introduced a colour cut. Given 
that a typical ~ 0.6 Mọ white dwarf with an effective temperature 
hotter than 5500 K has Ggp — Grp < 0.86 mag, we restricted our 
subsequent analysis to objects in this range of Gaia colour. 

Among the 2886 selected white dwarfs within 100 pc 
with spectral clasification in the MWDD, there are 2400 with 
Gpp — Grp < 0.86 mag, that we adopted as a validating sample, 
represents ~ 30% of our 100 pe white dwarf sample in this colour 
range. In particular, the MWDD validating sample contains 1789 
white dwarfs classified as DAs, and 611 white dwarfs of other spec- 
tral types DB, DC, DQ, DZ, and similar types, that we assigned to 
the non-DA group. Thus, the MWDD sample has a ratio of DA to 
non-DA of 74:26 with a significant fraction of the different spectral 
types. This fact, together with the aforementioned 30% cross-match 
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Table 3. Summary of the white dwarf spectral estimator built for the present study. 


Estimator name Reference Atmosphere models mean number points per SED 
VOSA This work Pure hydrogen (Koester 2010) 20 
Nearly pure helium H/He= -6 (Koester 2010) 
VOSA-GJP This work idem 56 
VOSA-GJP-trunc This work idem 56 
GF21-I Gentile Fusillo et al. (2021) Pure hydrogen (Tremblay et al. 2011; Kowalski et al. 2010) 3 
Pure helium (Bergeron et al. 2011) 
GF21-II Gentile Fusillo et al. (2021) Pure hydrogen (Tremblay et al. 2011; Kowalski et al. 2010) 3 
Mixed H/He=-5 (Tremblay et al. 2014; McCleery et al. 2020) 
MON22 Gaia Collaboration et al. (2022a) * 60 
VOSA VOSA-GJP VOSA-GJP-trunc 
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Figure 4. Probability distribution of being DA for the different estimators analyzed in this work for the MWDD validating sample. 


with our catalogue, places the MWDD as an excellent validating 
sample for our spectral estimator tests. 


4.3 Results: confusion matrices and scores 


In Fig. 4, we show the probability distribution of being DA for the 
set of estimators under study. The ordinate axis represents the num- 
ber of white dwarfs, while blue and red histograms correspond to 
white dwarfs belonging to the DA and non-DA classes, respectively, 
of the MWDD validating sample. A first glance to the distributions 
reveals that the pair of estimators VOSA-GJP and VOSA-GJP-trunc 
(top-middle and top-right panels, respectively) have a similar be- 
haviour, indicating an excellent capability for disentangling DA 
from non-DA objects. In the case of the VOSA and MON22 esti- 


mators (top-left and bottom-right panels, respectively), the general 
behaviour is also acceptable. However, the fraction of non-DA white 
dwarfs misclassified as DAs is larger than in the two previous esti- 
mators. The last pair of estimators, GF21-I and GF21-II (bottom-left 
and bottom-middle panels, respectively), seem to have a similar be- 
haviour, correctly identifying practically all DA white dwarfs, but 
failing in classifying most of the non-DA objects. 


The previous analysis can be quantified by means of the con- 
fusion matrix, which represents the correlation between the true 
label object (in our case the MWDD DA or non-DA label; rows) 
as a function of the predicted label (the DA or non-DA estimator 
label; columns). The obtained results are shown in Fig. 5, where 
the displayed values correspond to the total number of objects, the 
normalized values (in brackets) for each class, and the colour is 
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Figure 5. Confusion matrices for the different estimators analyzed in this work. Displayed values correspond to the total number of objects, while in brackets 
the normalized value with respect to the total population. Colour is scaled proportional to the total number of objects considered 


scaled proportional to the total number of objects considered. The 
ideal case would correspond to a diagonal matrix with a fraction of 
DA:non-DA of 74:26, that is the fraction presented in the MWDD 
validating sample. Our estimator VOSA-GJP actually reproduces 
this ideal case, showing an excellent performance: only 6% of true- 
DA white dwarfs are misclassified as non-DA, while less than 17% 
of true non-DAs are not properly identified. 

In Table 4 we list the accuracy and F1 global scores, as well 
as the sensitivity (also named recall) and the precision for each of 
the analyzed estimators. As previously stated, the estimator VOSA- 
GJP presents an excellent performance in all the analyzed scores. 
Particularly good are the accuracy and the Fl-score, 0.91 and 0.94 
respectively, as well as the precision of 0.94, reaching these scores 
the highest values among all the estimators. Regarding the sensi- 
tivity, the VOSA-GJP-trunc and GF21-II estimators, with 0.98 and 
0.95, respectively, present a slightly better performance than VOSA- 
GJP, 0.94. However, the better ability of these estimators to retrieve 
objects of a certain class is worsened by a substantial lower perfor- 
mance regarding the precision with which an object is identified. All 
in all, the VOSA-GJP can be considered as an excellent estimator 
for the DA and non-DA white dwarf identification. 

Finally, before analyzing the physical properties of the identi- 
fied samples in the next section, a few more checks have been made. 
As previously stated, for our initial analysis we adopted a thresh- 
old value of 7=0.5. Now, we left 7 as a free parameter. For the 
VOSA-GPJ estimator, the value that maximizes the performance is 
still 7 = 0.5. For the rest of estimators, slight variations around this 
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Table 4. Summary of the performance score, sorted from the highest to the 
lowest accuracy, for the different estimators under study. 


Estimator Accuracy Fl-Score Sensitivity Precision 
VOSA-GJP 0.91 0.94 0.94 0.94 
VOSA-GJP-trunc 0.88 0.92 0.98 0.88 
MON22 0.78 0.85 0.81 0.89 
VOSA 0.76 0.85 0.87 0.82 
GF21-II 0.70 0.82 0.96 0.72 
GF21-I 0.57 0.68 0.64 0.72 


value are found. However, in none of the cases the performance of 
the rest of estimators improved that of the VOSA-GPJ. 


5 THE PHYSICAL PROPERTIES OF DA AND NON-DA 
100 PC WHITE DWARF POPULATION 


In this section, we analyze the physical properties of the 100 pc 
white dwarf population once classified according to DA or non-DA 
spectral type by our estimator VOSA-GJP. 

We recall here that our analysis is restricted to those white 
dwarfs with good model fit to the GJP SED and in the Gaia colour 
range Gpp — Grp < 0.86 mag, which approximately corresponds to 
a typical ~ 0.6 Mọ white dwarf with an effective temperature just 
above 5500 K. Thus, we did not expect BP-band photometry to 
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Figure 6. Gaia Hertzsprung—Russell diagram of all Gaia sources consid- 
ered in this work. Those with an available GJP SED and a reliable model 
fit (Vgfp < 15) are marked in blue. These objects in the range of colour 
Gpp— Grp < 0.86 (dashed line) are thoroughly analyzed in Section 5. Those 
sources with a bad fit are shonw in green, and those with no Gaia spectrum 
available in red. The total percentage of excluded objects represents only 
1% in the considered colour range. 


be affected by any bias in this colour range, so we use this more 
commun Gaia colour in the analysis hereafter. 

We show in Fig. 6 the Gaia HR diagram for our 100-pc sample, 
where those objects with an available GJP SED obtained from the 
non-truncated Gaia spectra (see Sect. 3.3) and a reliable model fit 
(Vgfp < 15) are marked as blue dots. These objects which represent 
the 97% of the entire 100 pc white dwarf population. The remaining 
7% of objects have either a poor GJP SED fit (green dots) or even 
no Gaia spectrum is available (red dots). As it can be seen, the 
vast majority of them are faint and red objects. For the colour range 
considered in the analysis (left to the dashed line), the total number 
of objects classified is 8150, which represents the 99% of the total 
population for that range of temperatures. This is up-to-date the 
most complete spectral type classified sample of white dwarfs. 


5.1 The Hertzsprung-Russell diagram for DA and non-DA 
white dwarfs 


In Fig. 7, we show the Gaia HR diagram for our sample of DA (top 
panel) and non-DA (bottom panel) white dwarfs. For comparative 
purposes, we have also plotted different spectral types of already 
classified white dwarfs according to the MWDD. The A, B and Q 
branches are also pointed (see Gaia Collaboration et al. 2018). 

A first glance at Fig. 7 reveals an excellent agreement between 
the MWDD classified white dwarfs and those classified by our best 
estimator. This is not surprising since the accuracy of the VOSA- 
GJP estimator is larger than 90%. Secondly, the loci in the colour- 
magnitude diagram of the DA and non-DA populations are clearly 
distinct. The characteristic bifurcation on the HR diagram has been 
put into manifest since the Gaia-DR2. However, the association 
of each of the two branches, A and B, to a pure-hydrogen or pure- 
helium atmosphere models, respectively, is a much debated problem 
(e.g. Jiménez-Esteban et al. 2018; Bergeron et al. 2019). In partic- 
ular, if we assume that the B branch, within 0.1 < Ggp — Grp < 0.5 


mag, is formed by white dwarfs with pure-helium atmospheres, 
this would correspond to average masses larger than the canonical 
0.6 Mo (e.g. Bergeron et al. 2019). Several alternative explanations 
to avoid this problem have been proposed, such as mixed hydrogen- 
helium atmospheres, or spectral evolution from hydrogen to helium 
envelopes, among others (e.g. Bergeron et al. 2019; Ourique et al. 
2020). In any case, the characterization of the A and B branches 
has been limited to the spectroscopically identified white dwarfs in 
that region, leading to a large fraction of unidentified objects and, 
consequently, to ignore the real fraction of DA and non-DA in these 
branches. However, through our analysis, we can confirm that the 
B branch is mostly but not exclusively formed by non-DA white 
dwarfs, i.e. 65% of non-DAs and 35% of DAs, while the A branch 
is practically formed by DA white dwarfs (less than 6% are non- 
DAs). Given that the proportion of sources in the A and B branches 
is 62 to 38%, respectively, this implies that non-DA white dwarfs 
represent ~ 25% of the objects, considering both branches together. 

A closer look to the HR diagrams reveals that DA white dwarfs 
are mostly following the main cooling track (A branch), correspond- 
ing toa canonical 0.6 Mọ white dwarf. However, as has been already 
mentioned, 18% of the objects follows in parallel the main branch 
along the B branch. This bimodality was pointed out by El-Badry 
et al. (2018) and it was associated to a flattened in the initial-to- 
final mass relationship, which causes an overproduction of 0.8 Mo 
white dwarfs. However, other possibilities may contribute to this 
issue, such as the contribution of white dwarf mergers (e.g Kilic 
et al. 2018). Furthermore, a third branch, starting in the Q branch 
region, is growing in number of objects and joining the main branch 
at Gpp — Grp ~ 0.8 mag. As we will analyze in Section 5.2, the net 
effect is an increase of the average mass as objects are getting cooler. 
The combined effects of crystallization, 22Ne sedimentation, white 
dwarf mergers, and probably some other delaying physical cooling 
processes, may lead to the formation of this extended Q branch (e.g. 
Cheng et al. 2019; Tremblay et al. 2019; Kilic et al. 2020; Camisassa 
et al. 2021; Blouin et al. 2021). 

Regarding the non-DA distribution, apart from the commented 
fact that they mostly are located in the B branch, we observe a widen- 
ing of the track for the coolest objects, i.e., Ggp — Gpp > 0.6 mag. In 
this case, the reverse of the effect than in the DA sample is presented: 
a decrease in the average mass is expected for cooler non-DA ob- 
jects. This effect, extended for cooler and less massive objects, was 
discussed in Bergeron et al. (2019), although no final conclusion 
was reached. In any case, the true nature of non-DA objects that 
constitute the B branch, as well as the peculiarities in both diagrams 
already commented, are beyond the scope of the present work. 


5.2 Stellar parameter distributions 


Relying on updated white dwarf evolutionary models from La Plata 
group and atmosphere models from Koester (2010), we estimated 
stellar parameters for the white dwarfs with GJP SED, a reliable 
model fit (see Sect. 3.3), and Ggp — Grp < 0.86 mag. We used dif- 
ferent models depending on the spectral classification of the source 
(see Sect. 4). For DA white dwarfs, we used the cooling models of 
Althaus et al. (2013) for low-mass helium-core white dwarfs, the 
models of Camisassa et al. (2016) for average-mass carbon-oxygen- 
core white dwarfs, and the models of Camisassa et al. (2019) for 
ultra-massive oxygen-neon-core white dwarfs. For non-DA white 
dwarfs, we employed the hydrogen-deficient cooling models of 
Camisassa et al. (2017) for average-mass carbon-oxygen-core white 
dwarfs and Camisassa et al. (2019) for oxygen-neon ultra-massive 
hydrogen-deficient white dwarfs. All these models include realistic 


MNRAS 000, 1-17 () 


100 pc WDs classification from Gaia-DR3 and the VO 11 


9 
, Gaia 100 pc sample 
10 This work DA 
i MWDD DA 
he "KE a ey MWDD DAH 
11 RE. tf fe. MWDD DAZ 
12 


-0.4 -0.2 00 02 04 06 08 1.0 
Gap ~G 


BP “RP 


Gaia 100 pc sample 
This work non—-DA 
MWDD DB 

A MWDD DBA 
MWDD DC 

MWDD DQ 

MWDD DZ 


-0.4 -0.2 00 02 04 06 08 1.0 
GG 


Figure 7. Gaia HR diagram for the DA (top panel; blue dots) and non-DA (bottom panel; red dots) white dwarfs classified in this work. Also plotted, for 
comparative purposes, our entire Gaia 100-pc sample (grey dots) and those classified in several spectral types in the MWDD. The A, B and Q branches are 
pointed for helping in the discussion. 
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initial chemical profiles that are the result of the calculation of the 
progenitor evolution and consider all the relevant energy sources 
that govern their evolution. That is, they include neutrino loses, the 
gravothermal energy released by the ions, the energy released by 
slow 22Ne sedimentation and the energy released during the crystal- 
lization process, both as latent heat and due to the phase separation 
(see Camisassa et al. 2022, for details of its implementation). We 
used the atmosphere models of Koester (2010) to turn these evo- 
lutionary models into the magnitudes in the Gaia passbands. For 
DA and non-DA white dwarfs we employed pure hydrogen and pure 
helium atmospheres, respectively. Considering the Gaia magnitude 
G and colour Gpp — Grp of each of the white dwarfs in our sample, 
we have interpolated in the theoretical models to obtain their masses 
and effective temperatures. We excluded all DA white dwarfs with 
estimated masses below 0.239 Mo and all non-DA white dwarfs 
with masses below 0.51 Mo, to avoid extrapolation uncertainties. 

In Figure 8, we show the distribution of white dwarf masses 
versus effective temperature for our sample of DAs (left panel) and 
non-DAs (right panel). Also marked are the lines of crystalliza- 
tion onset assuming CO-core for all white dwarfs (solid lines) or 
adopting ONe-cores for masses above 1.05 Mọ (dotted lines). The 
crystallization onset was specifically calculated for DA (magenta 
lines) and non-DA (blue lines) white dwarfs. The distributions ob- 
tained reflect most of the peculiarities previously discussed regard- 
ing the Gaia HR diagram and already analyzed in similar diagrams 
by Bergeron et al. (2019) and Kilic et al. (2020). The bulk of the 
DAs is formed by objects of ~0.6Mọ uniformly distributed for 
the full range of temperatures. A secondary similar group, with 
masses around ~0.8 Mo but more diffuse can be intuited. More- 
over, following the DA crystallization onset line, a group of massive 
and hot white dwarfs (starting with masses from ~ 1.2-1.3 Mọ and 
~ 30 000 K), whose average mass is decreasing with the temperature 
is clearly visible. This last trend of stars corresponds to the Q branch 
and its extension, and they are related to the slow down on the white 
dwarf cooling rate due to the energy released by crystallization and 
by the sedimentation of 2”Ne (Tremblay et al. 2019; Camisassa et al. 
2021). The combined effect of these two later groups of stars for 
temperatures below 9000 K is an increase on the canonical 0.6 Mo 
average mass. Lastly, an increasing number of low-mass stars (due 
to double-degenerate or unresolved binaries) for lower temperatures 
is also observable. 

Regarding the non-DA distribution of masses and tempera- 
tures, we found a lack of ~0.6 Mo objects for temperatures be- 
tween 7000 < Tag < 10000 K. This effect, as pointed out in Berg- 
eron et al. (2019), can be avoided with the use of other atmospheric 
models such as mixed hydrogen-helium envelopes, instead of the 
pure-helium models we used. Likewise, the bulk of white dwarfs 
that have an average mass of ~0.6Mọ at ~ 20000K, presents a 
smaller value of ~ 0.55 Mo at lower temperatures of ~ 6000 K. In 
any case, it is out of the scope of the present study to ascertain the 
ultimate atmospheric composition of these objects. Independently 
of this issue, the non-DA distribution also shows a clear increase in 
the number of massive objects for cooler temperatures. 

An individual analysis of the parameters is carried on in Fig. 9, 
where we show the frequency distribution of effective temperature 
(left panel) and mass (right panel) for our samples of identified 
DA (blue histogram) and non-DA (red histogram) white dwarfs. 
Roughly speaking, the distributions of effective temperatures seem 
to be similar for DA and non-DA white dwarfs. Although a deficit 
of non-DA objects is observable for hotter stars, Tag Z 12 000 K, 
an excess occurs for cooler temperatures, Tag £9000 K. We further 
investigate this issue in the next section. 


Regarding the mass distribution, the discrepancies are notori- 
ous. The DA distribution shows the canonical peak at = 0.58 Mo, 
a bump extended up to 0.8 Mo, and a small fraction of low-mass 
(non-single origin) white dwarfs. These trends are in agreement 
with previous reported works (e.g. Jiménez-Esteban et al. 2018; 
Kilic et al. 2020). With regard to the non-DA distribution, a clear 
bimodality with peaks at ~0.6 Me and ~ 0.76 Mo is present. As 
stated, this second peak is a consequence of the bifurcation in the 
HR diagram and the use of pure-helium atmospheres to derive its 
mass. Otherwise, for instance, in the case of mixed hydrogen-helium 
atmospheres, this second peak can eventually be removed, or at least 
reduced. However, if we adopt a flattened initial-to-final mass rela- 
tionship as in the case of the DA population, a similar peak around 
0.8 Mo should be expected in the non-DA mass distribution. Under 
this hypothesis, and assuming the same ratio between the peaks at 
0.6 and 0.8 Mo as in the DA case (i.e., 80:20 respectively, see Fig. 9 
right panel), we would expect that 20% of the non-DA white dwarfs 
in the B branch to have masses in the range of ~ 0.8 Mo. 

Further information can be retrieved when the mass distribu- 
tion is depicted as a function of the effective temperature. This 
is done in Fig. 10, where we plot the DA (left panels) and non- 
DA (right panels) mass distribution for hot (Tag >12000 K; top 
panels), medium (7000 < Tag < 12 000 K; middle panels), and cool 
(Te < 7000 K; bottom panels) white dwarfs. These temperature 
ranges approximately correspond, respectively, to the regions be- 
fore, during, and after the bifurcation in the Gaia HR diagram. 

The DA mass distribution for hotter white dwarfs presents the 
canonical peak at = 0.58 Mo and an extended tail for massive white 
dwarfs. However, for Tag <12000K practically no objects more 
massive than 1.0 Mo are found. Two factors should be taken into 
account. On one hand, as it is seen in Fig.6, the 1% of sources 
excluded (red and green dots) in the analysed colour range have 
larger absolute magnitudes and redder colour. These objects typi- 
cally correspond to cool massive white dwarfs. On the other hand, 
as pointed out by Kilic et al. (2020), this absence of more massive 
stars can be also associated to the delay in the cooling process due 
to the effect of core crystallization, among others. Moreover, while 
the main peak remains almost constant regardless of the tempera- 
ture range, the second one at ~ 0.8 Mo seems to be more prominent 
in the middle range of temperatures. Even more, for the coolest 
white dwarfs an extended tail from 0.6 to 0.8 Mg appears. This 
feature is indicative that, even assuming the bimodality is caused 
by the flattened of the initial-final mass relationship, an additional 
process such as the contribution of the merger white dwarf popu- 
lation needs to be added to this scenario. Some of these features 
on the mass-temperature distribution have been partially guessed 
in previous works (e.g. Kepler et al. 2007; Tremblay et al. 2010; 
Rebassa-Mansergas et al. 2015a,b; Tremblay et al. 2016). However, 
it was not until the arrival of the data provided by the Gaia-DR3 
that we had a full picture of the mass-temperature distribution of 
the white dwarf population. 

In the case of the non-DA mass distribution, the dependence 
with the effective temperature is stronger than in the previous case. 
While the hottest non-DA white dwarfs peak at = 0.58 Mọ with an 
extended tail up to massive white dwarfs, the non-DA mass distri- 
bution for medium temperatures peaks at = 0.75 Mo. We recall that 
this range of effective temperatures corresponds to the bifurcation 
and that the increase of the average mass is likely a consequence 
of the use of pure-helium atmospheres for white dwarfs in the B 
branch. In the case of the coolest non-DA white dwarfs, an extended 
fraction of low mass (< 0.6 Mọ) white dwarfs appears. Neverthe- 
less, all these features on the non-DA mass distribution should be 
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Figure 8. White dwarf masses vs effective temperature for our sample of DA (left panel) and non-DA (right panel). In addition, the crystallization onset 
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Figure 9. Distributions of effective temperature (left panel) and mass (right panel) for our sample of identified DA (blue histograms) and non-DA (red 


histogram). See text for details. 


taken with caution, given that we are assuming pure-helium atmo- 
spheres for all non-DA objects, and that a combination of different 
atmospheres such a mixed, carbon-contaminated, etc., is expected 
to be present in the observed sample. 


5.3 The ratio of DA to non-DA white dwarfs 


The ratio of DA to non-DA white dwarfs as a function of the effective 
temperature, also commonly referred to as the spectral evolution, 
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is a capital tool for understanding the evolutionary processes in 
white dwarf atmospheres (e.g. Ourique et al. 2020; Cunningham 
et al. 2020, and references therein). Due to the diversity of selec- 
tion effects, one of the major problems that previous studies had to 
face when building a detailed spectral evolution function was the 
incompleteness of the analyzed samples. However, based on the ro- 
bustness of our estimator built from Gaia spectra, we could analyze 
with a highly degree of completeness our statistically significant 
100 pc white dwarf sample. That is, we can consider it as a volume 
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Figure 10. Distribution of mass for our sample of identified DA (left panels) and non-DA (right panels) in different range of temperatures. 


limited sample, thus avoiding the need of any assumption on cor- 
recting factors. Furthermore, our analysis allowed to reconstruct the 
spectral evolution for a wide range of effective temperatures, from 
25 000 to 5500 K. 


In Fig. 11, we show the ratio of non-DA white dwarfs with 
respect to the total number of objects as a function of effective 
temperature (black line). We recall that our classification was done 
for the Gaia colour range Ggp — Grp <0.86 mag, corresponding 
to Te >5500K and, in order to avoid the effect of unresolved 
binaries, we selected objects only with Mwp > 0.51 Mo. Error bars 
were estimated as of = yf x (1 — f)/N, where f is the fraction 
of non-DA to the total number of objects, N, and weighting the 
contribution of each object according to the probability of being DA. 
For comparative purposes, we also show some other distributions 
recently published in the literature. 


The first characteristic of our spectral evolution function is the 
general smooth behaviour, a direct consequence of the high degree 
of completeness of our sample. Secondly, error bars are larger for 
hotter temperatures, given that our 100 pc volume-limited sample is 
relatively small in comparison to magnitude-limited samples, and 
thus, the scarcity in hot objects. However, the higher resolution of 
our distribution for cooler temperatures is evident. Analyzing the 
particular details of our distribution, we discern some of the char- 
acteristics already presented in previous works. For instance, in the 
range of temperatures between ~ 23 000 and ~ 13 000K, our spec- 
tral distribution is nearly constant with an average ratio of ~ 13% and 
a slight growth up to 13 000K. This behaviour is compatible with 
most of the distributions: Ourique et al. (2020), Lopez-Sanjuan et al. 
(2022), Cunningham et al. (2020), and Genest-Beaulieu & Bergeron 
(2019). At that point, there is an abrupt increase that reaches its 
maximum fraction, f = 36% at Teg ~ 10500 K, which corresponds 


to the bifurcation zone in the Gaia HR diagram. A similar high 
ratio is also found in Cunningham et al. (2020); Ourique et al. 
(2020). In particular, our distribution perfectly resembles that of 
Cunningham et al. (2020). For cooler temperatures, the fraction 
decreases down to f £0.31 for the coolest bin (Tag ~ 5500 K). A 
similar trend is also found by McCleery et al. (2020). On the con- 
trary, we have found no evidences of an increase ratio of non-DA 
to DA between 6250 < Top < 7500 K due to convective mixing, as 
predicted by Blouin et al. (2019). However, as it is evident from 
the discrepancies among the distributions, any conclusion in the 
cooler range of temperatures should be taken with caution, as larger 
photometric errors and selection biases dominate in this region. 


6 CONCLUSIONS 


An update of our 100 pc white dwarf sample (Jiménez-Esteban 
et al. 2018) was performed, taking advantage of the new photomet- 
ric and astrometric data together with the low resolution spectra 
provided by the Gaia-DR3. A total of 12718 white dwarfs config- 
ure our new catalogue, from which 12 342 have Gaia spectra. This 
constitutes the largest nearly-complete volume-limited white dwarf 
sample available to date. 

The use of automated algorithms for fitting and analysing 
SEDs, provided by the Spanish Virtual Observatory, and in par- 
ticular VOSA, allowed us to extract the maximum information from 
the SED of the white dwarfs in our sample. In this sense, we built 
an estimator of the probability of being DA from the fitting of 
hydrogen-rich and nearly helium-pure atmosphere models to the 
synthetic J-PAS photometry derived from the Gaia spectra, when 
all the coefficients are taken into account. The statistical analysis 
carried out when comparing this estimator with two others also 
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Figure 11. Gaia Ratio of DA vs non-DA for our sample of identified white dwarfs (black points and lines). Objects with Mwp < 0.51 Mọ has been disregarded 
to avoid contamination from unresolved binaries. For comparative purposes we also show spectral evolution distributions from other works. 


built in this work and others available from the literature, revealed 
its superior performance. The validating test by using the spectral 
classification from the Montreal White Dwarf Database showed an 
accuracy larger than 90%, together with a precision and a sensitivity 
of 94%. 

The excellent results in the identification of DA and non-DA 
white dwarfs led us to apply this estimator to the entire population 
with an effective temperature above 5500 K. A total of 8150 objects 
(representing the 99% of objects of our catalogue in that range 
of temperatures) have been spectrally classified in these two main 
groups. This classification has allowed to precisely determine the 
proportion of DA and non-DA objects in the different regions of the 
Gaia HR diagram. 

Our results showed that the A branch in the region 0.1 < Ggp — 
Grp < 0.5 mag is practically formed by DA stars (94%), while the 
B branch, although mainly constituted by non-DA white dwarfs 
(65%), contains a significant fraction of DA white dwarfs (35%). 
These results imply, on one hand, the confirmation of a bimodality 
in the DA mass distribution. Some hypothesis have been outlined 
such a stepper initial-to-final mass relation, or the contribution of 
white dwarf mergers (e.g. El-Badry et al. 2018; Kilic et al. 2018). 
On the other hand, the practically nonexistence of non-DA white 
dwarfs in the A branch is indicative that helium-pure atmospheres 
suffer from some mixing process, leading to mixed envelopes for 
temperatures below 12000 K. However, we cannot ensure that all 
non-DA white dwarfs in the B branch present a mixed atmosphere. 
Assuming the bimodality argument of the DA population caused 
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by the flattened initial-to-final mass relation, in the case of non-DA 
white dwarfs, this scenario would imply that 20% of non-DAs are 
genuine ~ 0.8 Mo white dwarfs. In any case, the ultimate charac- 
terization of these atmospheres is beyond the scope of the present 
work. 

Our analysis also allowed us to derive stellar parameters by 
means of photometric interpolation in updated cooling sequences 
from La Plata models with hydrogen-pure or helium-pure atmo- 
sphere models. This way we built, for instance, the mass distribution 
for DA and non-DA stars. Our results showed a bimodality in the DA 
mass distribution, as previously commented and already reported 
in the literature, with a main peak at ~ 0.58 Mo and a secondary 
one at ~ 0.8 Mo. An extended tail for white dwarfs more massive 
than 0.8 Mọ is also observed for Tog > 12000 K, that disappears 
for lower temperatures due to a combination of crystallization and 
other physical process along with incompleteness bias effects. It 
was also noticed that for temperatures below 7000 K the second 
peak fades. A possible cause is that, for this low range of tempera- 
tures, crystallization starts to happen for white dwarf less massive 
than 0.8 Mo, smoothing out the mass distribution between the two 
peaks. Regarding the non-DA mass distribution, although a prelim- 
inary guess was done by using helium-pure models, some conclu- 
sions may be drawn. A clear peak at ~ 0.58 Me mainly formed by 
genuine DB stars is present for temperatures above 12000 K. This 
peak is moved to ~ 0.72 M@ at lower temperatures as a consequence 
of fitting the B branch with helium pure atmospheric models. How- 
ever, for the lowest temperatures, the initial peak at even slightly 
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less massive value ~ 0.56 Mọ reappears. These facts indicate that 
individual spectral identification of white dwarfs atmosphere — i.e., 
DB, DBA, DQ, DZ, etc — is required in other to properly determine 
the non-DA mass distribution. 

Finally, our analysis also allowed us to derive a precise ratio of 
non-DA to DA white dwarfs as a function of the effective tempera- 
ture. Our spectral evolution distribution revealed a ~ 13% non-DA 
fraction nearly-constant smoothly increasing up to effective temper- 
atures of ~ 13 000 K, in excellent agreement with most of the recent 
published distributions. For lower temperatures, a marked growth 
was found reaching its maximum fraction of 36% at ~ 10 500K, in 
perfect agreement with the spectral evolution function of Cunning- 
ham et al. (2020) and also compatible with that from Ourique et al. 
(2020). Although for lower temperatures our result should be taken 
with caution, we found a decrease in the ratio of non-DA leading to 
a fraction of ~ 31% at ~ 5500 K. 

In summary, the conjunction of accurate white dwarf models, 
the excellent observational database provided by Gaia, and the ca- 
pabilities of automated tools such as the Virtual Observatory SED 
Analyser, allowed us to achieve the up-to-date largest spectral char- 
acterisation of the volume-limited white dwarf population in our 
solar neighbourhood, open thus the door to future studies such as its 
luminosity function or the star history of this region of the Galaxy. 


7 DATA AVAILABILITY 


The 100 pc White Dwarf Gaia-DR3 Catalogue is available online 
as supplementary material hosted by the jornal and at The SVO 
archive of Gaia white dwarfs? (see Appendix A) at the Spanish 
Virtual Observatory portal!®. 
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APPENDIX A: ONLINE CATALOGUE SERVICE 


In order to help the astronomical community on using our cata- 
logue of white dwarfs, we developed a wed archive system that can 
be accessed from a webpage!! or through a Virtual Observatory 
ConeSearch!?. 

The archive system implements a very simple search interface 
that permits queries by coordinates and radius as well as by other 
parameters of interest. The user can also select the maximum num- 
ber of sources (with values from 10 to unlimited) and the number 
of columns to return (minimum, default, or maximum verbosity). 

The result of the query is a HTML table with all the sources 
found in the archive fulfilling the search criteria. The result can also 
be downloaded as a VOTable or a CSV file. Detailed information on 
the output fields can be obtained placing the mouse over the question 
mark (“?") located close to the name of the column. The archive also 
implements the SAMP! (Simple Application Messaging) Virtual 
Observatory protocol. SAMP allows Virtual Observatory applica- 
tions to communicate with each other in a seamless and transparent 
manner for the user. This way, the results of a query can be easily 
transferred to other VO applications, such as, for instance, TOPCAT. 


This paper has been typeset from a TEX/IATEX file prepared by the author. 


ll http://svo2.cab.inta-csic.es/vocats/v2/wdw/ 


12 e.g. http://svo2.cab.inta-csic.es/vocats/v2/wdw/cs.php?RA=301.708&DEC=- 


67.482&SR=0.1&VERB=2 
13 http: //www.ivoa.net/documents/SAMP/ 


